Interacting fermions in two dimensions: beyond the perturbation theory 
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We consider a system of 2D fermions with short-range interaction. A straightforward perturbation 
theory is shown to be ill-defined even for an infinitesimally weak interaction, as the perturbative 
series for the self-energy diverges near the mass shell. We show that the divergences result from 
the interaction of fermions with the zero-sound collective mode. By re-summing the most divergent 
diagrams, we obtain a closed form of the self-energy near the mass shell. The spectral function 
exhibits a threshold feature at the onset of the emission of the zero-sound waves. We also show that 
the interaction with the zero sound does not affect a non-analytic, T 2 -part of the specific heat. 

PACS numbers: 71.10.Ay,71.10.Pm 



The Landau Fermi-liquid (FL) theory states that the 
low-energy properties of an interacting fermion system 
are similar to those of an ideal Fermi gas 0. A neces- 
sary (but not sufficient) condition for the validity of this 
theory is that various characteristic properties (thermo- 
dynamic parameters, quasi-particlc lifetime, etc.) can be 
expressed via regular perturbative expansions in the in- 
teraction. Back in the 50s, it was shown to be the case for 
a model of 3D fermions with both short- and long-range 
(Coulomb) repulsion. Soon thereafter, it was reazlied 
that the perturbation theory is singular in ID and, as a 
result, the FL is destroyed. The case of D = 2 had been 
the subject of an active and fairly recent discussion, with 
a number of proposals, most notably by Anderson |2j, for 
the breakdown of the FL in 2D. Although the prevailing 
opinion currently is that the FL is stable in 2D for suf- 
ficiently short-range (including Coulomb) interactions, a 
full description of the FL in 2D is still lacking. 

One of the problems in describing a 2D interacting sys- 
tem is that, in a contrast to the 3D case, a naive pertur- 
bative expansion in the interaction is singular. Namely, 
for a linearized single-particle spectrum, the imaginary 
part of the self-energy diverges on the mass shell. The 
divergence is logarithmic to second order in the interac- 
tion (3, S IE @ but, as we will show in this paper, it is 
amplified to a power-law, starting from the third order. 
This mass-shell singularity in 2D is a weaker form of the 
"infrared catastrophe" in ID. There, an on-shell fermion 
can emit an infinite number of soft bosons-quanta of 
charge- and spin-density fluctuations 7] , which gives rise 
to a power-law divergence of the self-energy on the mass- 
shell. In 2D, this mechanism is weakened but not com- 
pletely eliminated (the "memory" about the ID infrared 
catastrophe is erased completely only for D > 2). At 
a first glance, the breakdown of the perturbation theory 
confirms the conjecture that the FL is destroyed in 2D 
Q. However, all it really means is that in order to ob- 
tain physically meaningful results the perturbation the- 



ory must be re-summed even for an arbitrarily weak in- 
teraction 0,0- 

In this Letter, we report the results of an asymptot- 
ically exact re-summation of the divergent perturbation 
theory for D = 2. In addition to obtaining a closed and 
divergence-free form of the self-energy to all orders in 
the interaction, this procedure also allows one to identify 
the nature of the singularities in a perturbation theory 
as originating from the interaction between the fermions 
and the zero-sound (ZS) collective mode. At any finite 
order of the perturbation theory, the collective mode co- 
incides with the upper edge of the particle-hole (PH) 
continuum. In 2D, this degeneracy is strong enough to 
generate the divergences in the self-energy. Once per- 
turbations are summed up to all orders, the zero-sound 
mode splits off from the continuum, and the power-law 
divergences disappear. The remaining log-divergences 
are eliminated by restoring the finite curvature of the 
spectrum near the Fermi surface [5J, |6( . Thus the FL sur- 
vives. However, the resulting self-energy contains sev- 
eral non-perturbative features that signal a deviation 
from the standard FL behavior. First, the imaginary 
part of a non-perturbative, ZS contribution to the self- 
energy, Ezs (w, k), is a non-monotonic function of the 
"distance" to the mass-shell, A = u> — efe, and exhibits 
an anomaly at the threshold for emission of ZS waves. 
This anomaly gives rise to a non-Lorentzian shape of the 
spectral function. Second, as the theory is regularized 
by the difference in the ZS and Fermi velocities, which 
is small for a weak interaction, Ezs is strongly enhanced 
near the mass shell. On the mass shell, ReEzs is of the 
same - C/ 2 o->|w|-order, as the perturbative contribution. 
It has been argued that the non- analytic, perturbative 
self-energy ReE cx J7 2 cj|cj| gives rise to a non-analytic, 
T 2 -tcrm in the specific heat C(T) HHH- Ther efore, it 
becomes an issue whether the non-perturbative part of 
the self-energy also contributes to the T 2 -term in C(T). 
We show that this is not the case, as the thermodynamic 
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FIG. 1: Maximally divergent self-energy diagrams to second 
(a) and third (b) orders in the interaction. 



potential is not affected by the mass-shell singularities. 

In what follows, we present the results and outline the 
main steps. Detailed calculations can be found in the 
extended version of this communication 01 • 

To understand the origin of the mass-shell singularity 
in arbitrary dimension D, consider the self-energy to sec- 
ond order in the interaction (first two diagrams in FigJQ, 
and focus on small-angle scattering of fermions with al- 
most parallel momenta (forward scattering), as shown in 
Fig. [2a,. 

The imaginary part of the retarded self-energy on the 
mass shell is given by (for u> > 0) 

ImS F (w = e k ) oc ul J dnj dQQ - 1 J dO D 

x5 (0 + v F Q cos0)ImII(fi,Q), (1) 

where uq = U (0) m/27r is the dimensionless coupling con- 
stant, n(J7, Q) is the polarization bubble, and dOi = 
d0 2 = d-d, d0 3 = sintfcM. ImII(fi,Q) is the 
probability amplitude to generate a PH pair of frequency 
CI and momentum Q. For D = 1, the PH continuum 
consists of a single line il = VfQ and ImII(0, Q) oc 
5(0 — VfQ). A product of two <5-functions in Q yields 
a non-intcgrablc singularity, which is the origin of an 
infrared catastrophe @. For D > 1, the PH contin- 
uum is a broad band specified by < VfQ, and the 
fcrmion spectral function is softened by angle-averaging. 
For D = 3, the resulting self-energy is finite on the 
mass shell. However, for D = 2 the angle-averaged 
fcrmion spectral function still has a square- root sin- 
gularity: J d95(fl — vfQ cos-d) oc (vfQ — 1 1^ | ) 1 / 2 f or 

1/2 

VfQ > whereas ImII(f2, Q) oc f2/ (vpQ — has 
another square-root singularity at the continuum bound- 
ary. Merging of the two square-root singularities results 
in a logarithmic divergence of ImS on the mass shell. 
Near the mass shell (|A| < \u>\), ImE oc In |A|. 



FIG. 2: Forward(a) and backscattering (b+c) scattering pro- 
cesses. 



Higher-order diagrams contain higher powers of di- 
vergent PH bubbles; these bubbles are either explicit, 
as in the ring diagrams of Fig. or generated by 
integrating out the fermionic energy and momenta. 
One can verify that the square-root singularities ac- 
cumulate, and the higher-order diagrams for the self- 
energy diverge as powers of A -1 . At the n-th order, 
ImE F oc <w 2 (w/A)" /2_1 6I(cj/A), where 6(x) is the 
step-function. By Kramers-Kronig relation, these singu- 
larities generate similar power-law divergences in RcS^?. 
Hence, to find the self-energy in D = 2 one has to re-sum 
the perturbative series even for an infinitesimally small 
U. 

This summation can be carried out explicitly for small 
U by collecting diagrams with the maximum number of 
polarization bubbles at each order (to third order, such 
diagrams are shown in Fig.^). The result is 



SF(p) = i J G(p-q) 

U(0) 3U{0) 



4t/(0) - 2U 2 (0)n{q) 



l-E/(0)n(g) l + U(0)IL(q) 



(2) 



where q = (Q, Vl) and J q ■ ■ ■ = T^ fim / d 2 q/ (2tt) 2 .... 
The last two terms in Eq. correspond to the inter- 
action in the charge and spin channels, respectively; the 
first two terms contain unaccounted parts of first- and 
second-order diagrams. All terms in Eq. J2J have per- 
turbative contributions from the PH continuum. In ad- 
dition, the charge part contains a non-perturbative con- 
tribution from the ZS pole in the charge-channel prop- 
agator. Near the pole, [1 - c7(0)n(q)]~ 1 oc ugQ 2 /(0 2 - 
c 2 Q 2 ), where c = vf{1 + Uq/2 + . . . ) is the ZS velocity. 
Notice that the quanta of ZS are not free bosons: the 
residue of the pole vanishes at Q — 0, as it is required by 
the translational invariance of the system. Substituting 
the ZS propagator into Eq. we obtain for the corre- 
sponding contribution to the self-energy 
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where A* = Uqw/2, and functions Fp>j are the real and 



3 



F,(x) 



1x 2 



within the perturbation theory the interaction with the 
continuum is indistinguishable from the interaction with 
the zero sound. The logarithmic divergence in the self- 
energy is cut off by a finite curvature of the spectrum on 
a scale A ~ cu 2 /W, where W ~ Ep is the bandwidth. 
The power-law divergences in the PH contribution are 
cut at the same scale A* as in the Ezs, this time because 
an increase of £/(0)ImII near the boundary of the PH 
continuum reduces the effective interaction for A < A*. 
However, contrary to the zero-sound contribution, Epn 
is smooth at A ~ A* as no Cherenkov-type condition is 
involved. Near the mass shell, Epn reduces to 



FIG. 3: Scaling functions Fi (x) and Fr (x) 
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F(x) = (1+ -x)Vl - x 



3 ^ln 1 -±^, (4) 



correspondingly. Plots of Frj are presented in Fig. [3] 
Notice that < 0) = and F R (x > 1) = 0. 

For large x, i.e., far away from the mass shell, Fi(x) oc 
1/y/x, and ImEzs oc UqUj 5 / 2 / '\fK. This behavior repro- 
duces the result of the perturbation theory. Expanding 
further in 1/x, one obtains higher powers of I/A. How- 
ever, we see from (gj that the emerging singularity at 
A = is actually cut off at x ~ 1, i.e., at A ~ A*. 
For small x, Fi(x) — 3tt 2 x 2 /8, and ImEzs oc A 2 for 
A -C A*. At A = 0, ImEzs vanishes. Vanishing of 
ImEzs on the mass shell is a Cherenkov-type effect: be- 
cause the zero-sound velocity c > Vp, an on-shell fermion 
cannot emit a ZS boson. The A 2 -behavior of ImEzs in- 
between A = and A = A* tracks an increase of the 
phase space available for the emission of ZS bosons. We 
emphasize that as A* oc Uq, the crossover at A ~ A* 
could have not been obtained within the perturbation 
theory. Near A = A*, the derivative of ImE(A) diverges: 
lmSj(x) oc 1/vA" — A*. The physical meaning of scale 
A* can be understood by noticing that a Cherenkov-type 
constraint restricts the frequencies of emitted bosons to 
the range < fl < (A/A*) u (for u > 0). This con- 
straint is relevant for A < A*; for larger A, emission of 
bosons with frequencies in the entire interval, allowed by 
the Pauli principle (0 < ft < ui) is possible. 

The real part of the self-energy varies slowly near the 
mass shell as Fr{x -C 1) = 1 + x. Right on the mass 
shell, ReEzs oc UqU> 2 . Notice that ReEzs is of order 
Uq, although it is obtained by summing up terms of or- 
der Uq and higher. This enhancement is due to a non- 
perturbative cut-off of the mass-shell divergences in the 
perturbation theory. 

The remaining contribution to the forward-scattering 
part of the self-energy comes from the PH continuum. 
It contains a logarithmic divergence to second order and 
also power-law divergences to third and higher orders, as 
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for \uo\ <C ui c ; 
for \ui\ » lo c , 



UqE f 



where uj c 

The sum Ezs + Epn is the total contribution to the self- 
energy from forward scattering. Another contribution to 
the self-energy comes from the processes in which the 
fermions move initially in almost opposite directions, and 
then either scatter by small angle (Fig. 12b) or backscatter 
(Fig.|2t). We refer to both these processes as "backscat- 
tering". The backscattering part of the self-energy, Eb, 
is regular on the mass shell and, for weak interaction, 
needs to be evaluated only to second order Taking 
the result for Eb from Ref. |6( and adding it up with 
the results for Ezs and Epn, we obtain for the total self- 
energy on the mass shell 



ReE (w) 



uj \u>\ 
8E F 
uj \u)\ 
8E F 



ImEM = f^m^ 
2tt Ep M 



( U - [ U + U 2k F - UQU 2kF \) 

U2k F {uq - u 2 k F ) ; (5a) 

(5b) 



i 2 1 

max 



where u 2 k F = mil (2k F ) /2n and « max = max{u , u 2 k F }- 
The scaling function $ (a;) approaches the constant values 
of 1(1/2) + u 2kF (u 2kF —uq)/{2uq), for x > oo and x — » 0, 
correspondingly. 

The first term in the first line of Eq. (|5a|l is the 
non-perturbative contribution from forward scattering, 
whereas the term in square brackets is the backscatter- 
ing contribution. We see that near the mass shell both 
contributions are of the same order, i.e., the interaction 
with the zero-sound collective mode modifies significantly 
the perturbative result for the self-energy. For contact in- 
teraction (uq — u 2 k F ), the zero-sound contribution can- 
cels out with the backscattering one, so that the net ReE 
vanishes at the mass shell. At the same time, the effect 
of the interaction with ZS on the on-shell ImE is rather 
benign: all one has is a smooth crossover function inter- 
polating between two limiting values of the prefactor in 
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FIG. 4: A log-plot of the spectral function A(ui, k) as a func- 
tion of x — 2(oi — ek)/u 2 uj for fixed uj. L u = 2 and 7 = 0.1 
A kink at x — 1 is a signature of the interaction with the 
zero-sound mode. 



a familiar w 2 In |cj|-dependence. This is the consequence 
of the fact that the non-perturbative, ZS contribution to 
ImE vanishes on the mass shell. 

At the same time, away from the mass shell, the 
full self-energy and the the spectral function A(cu, k) = 
(— l/7r)ImG (uj, k) are affected by the non-monotonic 
variation of ImEzg with A near a threshold for emis- 
sion of ZS bosons. We consider a setup, when u> is fixed 
and the spectral function is measured as a function of 
k. This is equivalent to varying x = A/ A* at fixed ui. 
Near the mass shell, the spectral function consists of two 
terms A(x) = A qp (x) + Azs(x), where 



A qp {x) 



TT 2 U 2 Ep X 2 + 7 2 



is a quasi-particle contribution and 



Azs = 



F I {x) 



ir 2 u 2 E F (x 2 +7 2 ) 5 



(6) 



(7) 



is a ZS contribution. We introduced L u = In (Ep/u 2 \oj\) 
and 7 = \ll>\L uj /(2ttEf), and assumed for simplicity that 
wo = U2k F = u. The ZS contribution Azs(x) has a sharp 
maximum at x = 1, where function Fi(x) in Fig. [31 has 
a peak. This maximum gives rise to a kink in total A(x) 
(see Fig.0J. A similar consideration shows that the kink 
is present also for the Coulomb interaction. This kink 
should be detectable in photoemission experiments on 
layered materials and in momentum-conserving tunnel- 
ing between parallel layers of 2D electron gases ^3] • 

Finally, we consider an issue whether non-perturbative 
effects, discussed above, are relevant for the specific heat 
of two-dimensional fermions. We argue that they are 
not. The simplest way to see this is to use the def- 
inition, C(T) = -Td 2 ~/dT 2 , and expand the ther- 
modynamic potential 5 in powers of U. The thermo- 
dynamic potential can be evaluated in the Matsubara 



formalism, 
-(m/27r)[l 



in which the polarization operator Tl ri 



fim|/ V {v F Q) 2 + is regular at any Q m 
and Q. As a result, the Matsubara series for H converges 
for a weak, short-range interaction, and there is no need 
for re-summation of the perturbation theory for S. Eval- 
uating 3 to second-order in U yields 



SC (T) 



2 , 3mC(3) T 2 

hk F - u oU2k F ) — - — -g-- (8) 



This result coincides with C(T) obtained in Ref. [6| by 
finding the perturbative, backscattering part of the self- 
energy first and then using the relation between C(T) and 
E pj. The agreement between the two results shows that 
non-parturbative, forward-scattering self-energy does not 
affect C(T). 

For completeness, we also verified explicitly that there 
is no contribution to C(T) from forward scattering. 
To this end, one can use the relation between C(T) 
and the Green's function which, to first order in 
(lu — £fc) _1 E(o;, k), reads 



sc ^ = —Wr 

7T 01 



i / de k [ duju^- (9) 
T J_ OQ J„ oc dw 



xi5(u>- c k )ReE R (u), k) - -V — - — ImS fl (w, k) 

7T uj - e fc 



Here no(u>) is the Fermi function. It is crucial that both 
KeY, R (uj,k) and ItoS r (uj, k) are present in Eq. ©. The 
forward scattering self energy, Ep — Ezs + Eph, depends 
on k due to the presence of a non-perturbative scale A*, 
and therefore ImEp yields a finite contribution to C(T). 
Substituting the results for Ezs and Eph into Eq. J3J, 
we find that the forward scattering contributions to the 
specific heat from ReEp and ImEp cancel each other; 
hence, there is no contribution to C(T) from forward 
scattering. 
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